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We study by Quantum Monte Carlo simulations the phase diagram of lattice hard core bosons with 
nearest-neighbour repulsive interactions, in the presence of a super- lattice of adsorption sites. For a 
moderate adsorption strength, the system forms crystal phases registered with the adsorption lattice; 
a "supersolid" phase exists, on both the vacancy and interstitial sides, whereas at commensuration 
the superfluid density vanishes. The possible relevance of these results to experiments on ^He films 
adsorbed on graphite is discussed. 



I. INTRODUCTION 

Films of ^He adsorbed on strongly attractive substrates 
such as graphite, have been the subject of intense exper- 
imental and theoretical investigation, motivated by the 
remarkable variety of phases that such films display.^ ^ 
Although this research subject is relatively old, it has 
recently enjoyed a resurgence of interest in the context 
of a search for a supersolid phase of helium, namely one 
displaying simultaneously crystalline order and friction- 
less flow. That such a phase might exist in the second 
layer of an adsorbed ^He film on graphite, was first sug- 
gested by Crowell and Reppy over a decade ago but this 
contention has been recently brought back to fore.^^ In 
particular, it is suggested that a supersolid phase may oc- 
cur in the vicinity (or in correspondence) of a crystalline 
phase of the second adsorbed helium layer, registered with. 
the underlying substrate. The most recent, first princi- 
ples numerical studies of helium films on graphite, have 
yielded no evidence of such a supersolid phase, as no such 
registered crystal is observed.^ Nonetheless, the general 
issue of interplay between boson localization, possibly 
induced by an external pinning potential, and superflu- 
idity, remains one of general interest in condensed matter 
and quantum many-body physics. Moreover, theoretical 
predictions may also soon be tested experimentally, pos- 
sibly in a more controlled fashion, with ultracold atoms 
in optical lattices.^ 

That an external potential can significantly alter the 
phase diagram of a system of interacting bosons, giving 
rise to additional phases, is well known. For example, we 
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Figure 1: (Color online). Structure of the pinning potential on 
the triangular (left) and square (right) lattices. Filled circles 
represent lattice sites at which the pinning potential is worth 
— e. 



have recently shown how a disordering potential can give 
rise to a glassy phase, as well as induce superfluidity in 
systems that do not display it, in the absence of disorder.!^ 
In this work, we investigate theoretically, by means of 
Monte Carlo simulations, the possible existence of su- 
persolid phases of many-boson systems in the vicinity 
of crystalline phases stabilized by external periodic po- 
tentials (such as that due to the adsorption sites of a 
corrugated substrate), i.e., not present in the phase dia- 
gram of the system in the absence of an external poten- 
tial. Our study is based on a lattice model of interacting 
bosons, similar to that which has been utilized in previ- 
ous theoretical works, as a minimal model of the very 
nearly t wo-d imensional (2D) first few ^He adlayers on 
graphite.l^EIl 

Specifically, we consider here the (hard core) limit of 
infinite on-site repulsion, i.e., no more than one boson 
per site, and include a repulsive interaction (of strength 
V) between particles occupying nearest-neighbouring lat- 
tice sites. For sufficiently strong nearest-neighbour re- 
pulsion, the ground state of the system is known to be a 
crystal, at particle density p=l/2 ("checkerboard" solid) 
on a square lattice,^ and p=l/3 (2/3) on a triangu- 
lar lattice.!^ We include here a sublattice of attractive 
sites as well, acting as a strong pinning potential (see 
Fig. [T]). The sublattice is purposefully chosen not to 
correspond to any crystal structure which the system 
forms in the absence of an external potential. For suf- 
ficiently strong adsorption, additional crystalline phases 
can be expected to appear, registered with the adsorption 
sublattice; henceforth, we shall refer to these crystalline 
phases as commensurate^ the term incommensurate be- 
ing used instead to refer to the solid phases that the sys- 
tem forms in the absence of any external potential. For 
example, our choice of pinning potential on the triangu- 
lar lattice is such that particle density p=pc=^/^ corre- 
sponds to a commensurate crystal, while =1/3 to 
an incommensurate one. In other words, here the terms 
"commensurate" and "incommensurate" are with respect 
to the pinning potential. 

The purpose of this study is to provide a simple theo- 
retical framework to interpret experimental studies prob- 
ing for possible (commensurate) supersolid phases of he- 
lium films on graphite. Although we mostly discuss here 
numerical results obtained on a triangular lattice geom- 
etry, we have observed the same general physical be- 
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haviour on the square lattice as well. 

Our main findings is that supersolid phases exist on 
both the interstitial and on the vacancy side of a com- 
mensurate (registered) crystal. However, the superfluid 
density always vanishes as the density hits a value corre- 
sponding to either a commensurate or incommensurate 
crystal. In this sense, the pinning potential does not 
give rise to fundamentally new behaviour, with respect 
to what is observed in this model near and/or at incom- 
mensurate crystal phases, in the absence of any external 
potential. The vanishing of the superfluid response 
at crystal density, appears therefore to be a general hall- 
mark of any phase labelled as "supersolid" , occurring in 
a system of this type, i.e., in the presence of an external 
pinning potential. 

This paper is organized as follows. In Sec. [llj we de- 
scribe the model Hamiltonian. In Sec. Ill, we briefly 
review the methodology adopted in this work, while our 
numerical results are presented in Sec. [IV| Finally, we 
outline our conclusions in Sec. IVl 
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Figure 2: (Color online). Density p versus chemical potential 
/i for V=6 and three different pinning potential strengths, 
namely e=l (filled squares), e=2.5 (stars) and e=5 (filled cir- 
cles). Statistical errors are smaller than symbol sizes. Results 
shown are for a triangular lattice with L=24. 



III. METHODOLOGY 



II. MODEL 



The 2D hard core Bose model with nearest-neighbour 
interactions is expressed as follows: 

H = -t ^^{aldj + h.c.) + ^ fiifij - ^ /i^n^ . (1) 



We consider here a lattice (either triangular or square) 
of N = L X L sites, with periodic boundary conditions. 
The sums (ij) run over all pairs of nearest-neighbouring 
lattice sites, at (a^) is the Bose creation (annihilation) 
operator for a particle at site i, fii = aja^ is the local 
density operator, t is the hopping amplitude, V is the 
strength of nearest-neighbour repulsion, while jUi = ja — ei 
is a site-dependent chemical potential. Here, = e if the 
site belongs to the pinning sublattice (see Fig. [T]), zero 
otherwise, e > being the strength of the pinning poten- 
tial. As mentioned above, a hard-core on-site repulsion 
is assumed, limiting the occupation of every site to no 
more than one particle. 

The pinning potentials have been chosen for deflnite- 
ness to correspond to commensurate density pc'=l/4 for 
the triangular lattice, pc'=l/3 for the square lattice (Fig. 
[T]). No particular physical signiflcance should be ascribed 
to these choices, motivated only by the goal of making 
commensurate phases lower in density than the incom- 
mensurate ones, as would be the case for the second layer 
of ^He on graphite, if a commensurate crystal exists. It 
seems reasonable to expect that the basic physical con- 
clusions ought to remain unaffected by a different choice 
of pinning sublattice. 



We perform grand-canonical quantum Monte Carlo 
simulations to study the ground state properties of ([T]), 
using the Worm Algorithm in the lattice path-integral 
representation. As the details of this c ompu tational 
method are extensively described elsewhere,^^^^!! be- 
cause the calculations performed here are standard, we 
shall not review it here, and simply refer interested read- 
ers to the original references. 

The results shown here correspond to a temperature 
T sufficiently low (typically /3 = 1/T ^ L), so as to be 
regarded as essentially ground state estimates. In order 
to characterize the various phases, we compute the su- 
perfluid fraction p^, as well as the static structure factor: 
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where < ... > stands for thermal average. The presence 
of crystalline long-range order is signaled by a finite value 
of 5'(Q) for some specific wave vector, in the thermody- 
namic limit. For the triangular lattice, Q = (7r,27r/A/3) 
is a wave vector corresponding to the registered (com- 
mensurate) crystal at p=pc=l/4:, while Q = (47r/3,0) 
to an incommensurate crystal with p=p/=l/3, 2/3. On 
the square lattice, Q = (47r/3, — 27r/3) corresponds to a 
registered crystal at p=pc'=l/3, and Q = (7r,7r) to an 
incommensurate (checkerboard) crystal at p=p/=l/2. 

Unless otherwise specified, the results discussed below 
pertain to a triangular lattice geometry. We have carried 
out careful extrapolation of the results to the thermody- 
namic limit. In general, we observed that estimates ob- 
tained on a lattice of L x L=576 sites are identical, within 
statistical uncertainties, with the extrapolated ones. 
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IV. RESULTS 

We begin by discussing the existence of registered 
(commensurate) solid phases for sufficiently large val- 
ues of e, i.e., the strength of the adsorption potential. 
Fig. [2] displays the density computed as a function of the 
chemical potential for = 6, for three different pinning 
potential strengths, namely e=l (ffiled squares), e=2.5 
(stars) and e=5 (ffiled circles). The value of V is large 
enough for the incommensurate solid phases to exist, in 
the model without pinning potential. 

For a weak pinning potential, plateaus in the density 
appear only at p = p/ = 1/3 and p=2/3, i.e., in cor- 
respondence of the incommensurate phases. A disconti- 
nuity of the curve signals a ffist-order phase transition 
between a superfluid and the incommensurate crystal. 
Analogously to what observed in the model without pin- 
ning potential,^ for sufficiently large V the p(/i) curve is 
continuous on the interstitial side (p > p/), as a super- 
solid phase exists. 

As the strength of the pinning potential is increased, 
two additional crystalline phases appear, one at p = pc = 
1/4, the other at p=5/8. Henceforth, we shall focus our 
attention on the phase of density pc^ which is registered 
with the pinning potential. The other phase arises from 
the competition between the pinning potential and the 
nearest-neighbour repulsion, and the basic physics at or 
near this density is the same as near pc- As shown in Fig. 
[2j the p{p) curve displays no discontinuities on either the 
vacancy or the interstitial side of the commensurate crys- 
tal. This is evidence of vacancy- and interstitial-doped 
supersolids, as we discuss below. 

On performing a sufficient number of runs in the (V, e) 
plane, we have computed the phase boundary lines shown 
in Fig. [3j between a superfluid and crystal at the two den- 
sities pc and pj. The left part of the figure refers to pc- 
The system is superfluid for e < ec{V)^ where ec{V) is 
the minimum strength of the pinning potential for which 
a commensurate crystal is present, as a function of the 
strength of the nearest-neighbour repulsion V. For the 
commensurate phase ec{V) is monotonically decreasing 
with F, as the presence of a strong nearest-neighbour re- 
pulsion, which causes the appearance of the incommen- 
surate crystalline phase at p/ = 1/3 also favours the for- 
mation of a commensurate crystal at pc (in fact, ec{V) 
approaches zero as F ^ 00). On the other hand, the 
right part of Fig. |3] shows that the pinning potential sup- 
presses crystallization at density p/, i.e., a greater value 
of V is needed to stabilize the incommensurate crystal 
at p/ = 1/3 if the external pinning potential is present. 
This is due to the lattice mismatch of two competing 
crystalline phases. A sufficiently large value of e causes 
the incommensurate phase to disappear altogether. 

We now discuss the superfluid properties of the system 
near crystallization. We begin by examining the physics 
of the system near a commensurate solid phase. Fig. |4] 
shows the superfluid fraction ps (upper panel) and the 
static structure factor S'(Q = (tt, 27r/A/3)) as a function 
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Figure 3: (Color online). Ground state phase diagram of 
([1]). Symbols lie at computed phase boundaries between a 
superfluid (SF) and a commensurate crystal (CC) for p = 
1/4 (left), and between a superfluid and an incommensurate 
crystal (IC) at p = 1/3 (right). Statistical errors are smaller 
than symbol sizes. 
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Figure 4: (Color online). Superfluid density ps (upper panel) 
and static structure factor S(Q) (lower panel) in case of V = 
4.0, e = 5.0 and a wave vector Q = (tt , 27t / ^/S) . Statistical 
errors are smaller than symbol sizes. 



of the particle density. The choice of parameters, namely 
V=4 and e=5, corresponds to a situation in which the 
only crystalline phase that the system forms is the com- 
mensurate one, at a density pc- 

Both Ps and S{Q) are everywhere finite^ except at ex- 
actly pc where the superfluid response vanishes. The 
fact that p(p) is continuous everywhere, allows one to 
rule out coexistence of two phases (superfluid and crys- 
tal) possessing only one of the two types of order. Thus, 
based on its strict definition, one would have to conclude 
that this system is everywhere "supersolid" , except at 
commensurate density. However, such a denomination 
appears to be meaningful (if at all) only in the vicinity 
of Pc, where the physical character of the phase can be 
surmised to be that of a commensurate crystal doped 
with either vacancies or interstitials. Away from pc, the 
nature of the system is basically that of a fluid with a 
density modulation arising from the pinning potential.^^ 
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That at exactly p = pc the superfiuid fraction vanishes, 
is a significant result that warrants a few comments. Su- 
persolidity in model ([T]) on the triangular lattice (it is 
not present on the square lattice) , requires that a crystal 
be doped with interstitial particles, i.e., the superfiuid 
density of an undoped crystal is always zeio^ However, 
one might speculate that the lower density commensurate 
phase stabilized by the external potential might enjoy dif- 
ferent properties than the incommensurate one, which is 
the only one observed in the absence of a pinning po- 
tential. We find, however, that the superfiuid density 
always vanishes at pc^ on both lattice geometries consid- 
ered here. We have also repeated the same analysis for 
different choices of the parameters, including those for 
which both commensurate and incommensurate phases 
exists, but the presence of an incommensurate phase at 
higher density does not alter the physics of the system 
in the vicinity of pc, i.e., a "supersolid" phase exists on 
both the interstitial and vacancy sides, but not at com- 
mensuration. 

This is a result of potential experimental relevance, as 
studies of adsorbed ^He films on corrugated substrate 
on graphite, for which claims of possible supersolid be- 
haviour near commensurate density are made, can deter- 
mine the superfiuid response as a function of coverage. 
It appears from our results that, to the extent that ([T]) 
can be regarded as a reasonable qualitative model of a 
thin helium film on a corrugated substrate, the superfiuid 
signal must vanish at the coverage corresponding the oc- 
currence of a commensurate crystal, if one is to make a 
claim of a "supersolid" phase near or at commensuration. 

The physics of the system near the incommensurate 
crystal phase is the same as in the absence of an ex- 
ternal potential.1^ In particular, the superfiuid density 
again always vanishes at pj. Here too, one might have 
expected that the weakening of the incommensurate crys- 
tal caused by the pinning potential could give rise to a 
"softer" crystalline phase, capable of superflow. What is 
observed, however, is that as long as the incommensu- 
rate crystal exists, the superfiuid density at pi vanishes. 
There is always a first-order phase transition on the va- 
cancy side; on the other hand, on the interstitial side, 
depending on the value of V one may have a first-order 
phase transition to a superfiuid or a second-order phase 
transition to a supersolid. 

We now consider the regime in which both the com- 
mensurate and incommensurate phases exist, and explore 
the quantum phase transitions between the two, with the 
possible occurrence of intervening phases. There are two 
possible scenarios which have been observed: one is a di- 
rect transition from the commensurate to the incommen- 
surate crystal, through a first order transition. In other 
words, there is only a jump in the curve p{p) from pc to 
/?/, as shown in Fig. [5] This occurs roughly in a regime 
where e << V, and is a scenario that appears to apply 
to the first layer of helium on graphite,^ or to films of 
molecular hydrogen adsorbed on graphite, or other cor- 
rugated substrates. "'^'^^ 
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Figure 5: (Color online). Density p versus chemical poten- 
tial /i for the ground state of ([T]), with V = 8, e = 3. The 
density jump signals a first-order phase transition between a 
commensurate and an incommensurate crystal. 



The other scenario, occurring for e ~ is a second- 
order transition from the commensurate crystal at pc 
to a doped supersolid phase, and then, as the density 
is increased, to a superfiuid, followed by a first order 
phase transition from the superfiuid to the incommen- 
surate crystal at pj. The transition from supersolid to 
superfiuid is indicated by the change in slope of the p{p) 
curve. Again, we should note that this superfiuid phase 
with a density modulation arising from the pinning po- 
tential. 



V. CONCLUSIONS 

We have studied the ground state phase diagram of 
lattice hard core lattice bosons with a nearest-neighbour 
repulsion, in the presence of an external periodic poten- 
tial mimicking a superlattice of adsorption sites. We have 
utilized an exact numerical method. Our goal, in using 
such a simplified model, was obviously not that of achiev- 
ing a realistic description of an actual adsorbed film on 
a given substrate, but rather to gain qualitative under- 
standing of the different phases that one may be able to 
observe experimentally. We have carried out our stud- 
ies on triangular and square lattice geometries, as well 
as different choices of external potential periodicity, al- 
ways chosen to be not commensurate with the crystalline 
phase that the system forms in the absence of a poten- 
tial. The main findings are independent of the lattice 
geometry and/or potential periodicity. 

Our results show that in specific circumstances "su- 
persolid" phases (whether or not such terminology is ap- 
propriate, given the crucial role played by the external 
potential, is a matter of debate) can exist in the vicin- 
ity of commensurate crystals stabilized by the adsorp- 
tion potential. A distinctive signature of the occurrence 
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of such phases is the vanishing of the superfiuid density 
at commensuration. This seems to be a universal fea- 
ture of this type of system, one that we would expect to 
see in any experiment claiming observation of a super- 
solid phase of adsorbed films (e.g., of ^He) on substrates 
such as graphite. It is also worth mentioning that these 
predictions may also lend themselves to possible experi- 
mental verification by means of ultracold atoms in optical 
lattices.l^ 
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